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ABSTRACT 

This    report  examines   gas   dynamic   stabilization   of  high   pressure, 
highly  non-uniform  electrical    discharges   for  two  electrode   geometries. 
The   geometries   studied   are    the  multiple-point   to   plane   and    the  multiple- 
wire   to  plane.      A  successful    analytical   model    is    developed   utilizing   the 
anode  extension    concept,   which   describes    the    interactions   and   stabilizing 
effects   of   flow,    separating    the  effects   of   convection    and   turbulence,    in 
the   point-plane   geometry.      Comparisons   of   theory  and   data  are   presented, 
illustrating   actual    observed   and   calculated   breakdown    voltages   and 
currents    for  a  number  of   cases.      Flow  stabilized  high   power  density 
electric  discharges   are   applicable   to  high   power  electric  discharge    lasers, 
magnetohydrodynami c  power  generation,    and  other  new  systems   and  technolo- 
gies. 


VI  i 


I  .   INTRODUCTION 

A.   Background 

Electric  discharges  have  been  the  subject  of  intensive  study  for  many 
years.  Numerous  books,  papers,  reports,  etc.  have  been  published  which 
characterize  and  analyse  the  various  complex  and  interrelated  processes 
which  take  place  in  discharges  of  several  gas  compositions,  pressures,  and 
electrode  geometries.  The  primary  objective  in  most  of  these  studies  has 
been  to  attempt  to  relate  the  atomic  level  processes  to  the  observed 
macroscopic  discharge  properties  such  as  breakdown  voltage,  current  and 
overall  discharge  stability  for  use  in  engineering  analysis  and  application 
in  design.   Some  of  the  more  established  applications  of  this  type  of 
gaseous  electronics  include  glow  discharge  lighting,  the  high  voltage  arc, 
corona  discharge  analysis  around  power  lines  and  power  line  design,  high 
voltage  switching,  spark  plug  and  spark  gap  design,  thermal  printing 
devices  and  many  others. 

More  recent  applications  of  this  electric  discharge  technology  have 
been  generally  aimed  at  systems  which  are   geared  to  creating  very  high  power 
densities  within  the  discharge,  thrusting  physicists  and  engineers  headlong 
into  a  direct  confrontation  with  the  stability  problem,  not  just  with  the 
idea  of  "defining  the  envelope",  to  use  an  aeronautical  term,  but  to  develop 
ways  of  expanding  the  natural  discharge  stability  boundaries  through  any 
practical  means.  The  recent  high  power  applications  of  which  we  speak 
include  magnetohydrodynami  c  (MHD)  power  generation,  plasma  chemistry, 
electri c  discharge  convection  lasers,  E-beam  ionizer-sustainer  lasers  and 
others. 


As   a    result  of   this  growing    requirement    for   increasingly   higher   power 
densities,   many  analyses  have  been   done    in    the   stability  area   and  much 
insight  has   been   gained   in   the  many  modes   and   regimes    from  which   discharge 
instability  can  develop  to   cause  a   severe    local    breakdown    in   a   delicate 
balance  of  parameters.      One  of   the   successful    techniques  which   has    appeared 
in    the  engineering  of  high   power  density   systems    is   the   pulsed   system 
where-in    the  pulse  time    is   shorter   than    the  characteristic   times   of   the 
various    instability  modes  which  must   be  avoided.      The  E-beam   and   pulsed 
ioni zer-sustainer   lasers  are  examples   of   the  success   of   this    approach 
(Reference    1).      However,   present  E-beam  systems   are    large,   expensive   and 
delicate.      In   addition,    in   order  to  achieve   high   average   power  output 
from   the  pulsed  systems,  while   the   systems   can   be   scaled  up    in    size    to   any 
practical    limit,  once   the   size    limit  has   been    reached,    the  only  other 
alternative   is   to   increase   the  pulse    repetition   frequency    (PRF) .      But    this 
option    is   also    limited,    usually   by   a   characteristic  optical    cavity   clearing 
time  which   depends  on    cycling  out   by   flow   the  products   of  one   pulse  before 
initiating   the  next  one.      Failure   to  allow  sufficient    clearing    time  will 
result    in   cumulative   space   charge   and   gas   heating    instabilities  with 
resultant    low  power  breakdown. 

Due    to   the   inherent   limitations    cited  above  and   the    implicit    depen- 
dence of   those   systems  on   gasdynamics,    the  present    research   has   been   pur- 
sued  in  order  to  better  define   the   role  of   turbulence    in   a  high   pressure 
discharge  as   has   been  observed  by  many    researchers    (References   2,    3,    ^,    5 
and  6).      An   understanding  of   the   fundamental   modes  of  gasdynamic   coupling 
which  have  been   observed   to  affect  breakdown  voltage   and   current   may    lead 


to  new  ways  to  improve  present  high  power  density  systems  through  bulk  or 
electrode  boundary  layer  stabilization  techniques;  and  may  also  point  the 
way  for  development  of  sufficiently  high  power  self-sustained  CW  electric 
discharge  systems  which  are  inherently  more  compact,  rugged,  and  less  ex- 
pensive than  the  types  mentioned  above. 
B.        Present  Work 

Section  A  of  Part    II    (Fundamentals   of   Gasdynamic  Effects   on   a   Highly 
Non-Uniform  Electrical    Discharge)    of   this    report   develops    in    detail    a 
theoretical    model    for   the   calculation   of   breakdown    voltage   of   a  multi- 
point   to   plane    (points   positive)    electric  discharge   based  on    the   "anode 
extension"   concept    initially   proposed   in    Reference   7   and   further  discussed 
in    Reference    3.      The  model    incorporates    flow  and   turbulence   effects    and    is 
similar    in   physical    basis    to   the  "corona   cloud"  model    of  Hepworth,    Klene, 
and  Tozer  discussed    in    Reference   8   for  a   static    (non-flowing)    non-uniform 
electric  discharge.      The  present  model    successfully   separates    the  gasdy- 
namics    into   the   pure   convective  and  pure   turbulence  effects,    and    isolates 
specific  electrode   geometry  and  gas   parameters    in    individual    terms.      This 
facilitates    understanding  of   the   fundamental    relationships    between    the  gas 
properties,   electrode  orientation   and  geometry,    and    the    individual    gasdy- 
namic  coupling  mechanisms. 

This   section   also  presents   a  qualitative  discussion   of   the   breakdown 
characteristics  of   the   uniform  and  non-uniform  discharge   geometries   high- 
lighting  the  similarities    in   a   qualitative  and   semi-quantitative  manner. 
It    is    shown    that  Meek's   breakdown   criterion    for   the  non-flowing,    uniform 
field  and   the  equation   developed    in   Section  A    for   the  non-uniform,    flowing 
discharge   transform   into  very  nearly    identical    expressions    in    the  middle 
ground  of  a   flowing   uniform  field.      This    illustrates    the  general    nature  of 


the  equations  of  Section    A  and  provides   straightforward  physical    insight 
into   the  final    breakdown  mechanisms    in    the  geometric  extremes.      Radial 
diffusion   and   spreading    is   considered  as    related   to  the  mean   or  "drift" 
field   in   a   highly  non-uniform  geometry. 

Section   B   of   the   report  discusses   the  typical    discharge   section    used 
in   Section  A  with  the   three   rows  of  pins    replaced   by   three   thin  wires    (see 
Fig.    lb).      A  preliminary  analysis  which  was   done    in    the    terms   of   Reference 
9   basically    indicated   the  potential    for   improving   discharge    input   power  by 
a   factor  of  2-3  with    the  wire  anode.      Various   degrees   of  discharge    input 
power  enhancement   have  been   noted   in   similar  configuration   modifications 
in    connection  with  MHD   laser/power  generation   applications    (Reference    10). 

Recent  experimental    results    in    this   area   are   presented  and   discussed 
in    relation   to   the   points-plane   configuration.      To   date   due   to  power 
supply    limitations  no   improvement    in   performance  has   been   noted,   and  some 
recommendations    for  further    research    in    this   area   are   presented,   when   a 
larger   power   supply    is   available. 

Section   C  of  Part    II    discusses  electrode   design   geometry   and  electric 
field   distributions    for  optimum  electron   production    and   for   best    reduced 
electric   field   strength    (E/N  or   E/P)    for   laser   pumping. 

Part    I  I  I    presents   a   data   summary  of   current   and  past   point-plane  and 
wire-plane    laboratory  work,    reviews    the    theory   and  experimental    trends, 
indicating    important   correlations   and   features   of  the   anode   extension   model 


II.      FUNDAMENTALS   OF   GAS    DYNAMIC   EFFECTS   ON   A   HIGHLY 
NON-UNIFORM   ELECTRICAL   DISCHARGE 


A.        Prediction  of   Flow  and  Turbulence  Effects   on   Breakdown  Voltage 

1 .      Description  of   Discharge  Geometry   for  Analytical    Formulation 
A  significant   body  of   breakdown  data   has   been   generated    in   var- 
ious  highly  non-uniform  electric   field  geometries  with   positive   and   neg- 
ative polarities   at   the   highly   stressed  electrode.      The  point-to-plane 
geometry  with    the   points   positive,   however,    has   been   the   particular   case 
through  which   the   present   analytical   model    has   been    developed;    but    the 
model    is   sufficiently  general    that  the   basic  flow  and  electric   field 
interactions   described   should   also   be   useful    in    analysis,    interpreta- 
tion,  and   trouble   shooting  of  many  other  non-uniform  geometries    and 
systems  both  with  and  without   flow. 

The   specific  apparatus   used    in    the  current    research,    described 
in   detail    in    Reference    11    consists   of  multi-point    to   plane  electrode 
arrangement    (see  Fig.    la)    which    is  mounted    in    a  wind   tunnel    with   avail- 
able   test   section   velocities  on    the  order  of  200  m/s   of   ambient    pressure 
air.      Varying    degrees   of   turbulence  can   be   introduced  by   placing    turbu- 
lence generating   screens   upstream  of   the  electrode   set-up.      Two   flow 
modes   are   discussed   in    this    report,   namely, 

a.  Normal    Flow  -   directed   through    the   turbulence   screen,    past 
the  positive  multi-pin   anode,    through    the    interelectrode   gap,    and   thence 
through   the  grounded   cathode,   which    is   designed   to   approximate   a   plane 
with    low  blockage   to  permit  maximum  gas   flow   rates    (see   Fig.    2a). 

b.  Reverse   Flow  -   same  as   above  except    that    the  anode   and   cat- 
hode positions   are   reversed    (see  Fig.    2b)    and   the  turbulence   screen  was 
omi  tted. 
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While  our  earlier  effort  which    included   cross-flow   tests    (Fig. 
2c)    was   based  on    an   assortment  of   turbulence   screens  which   produced  a 
variety  of  spectra  and    intensities,    the   present  work    (which    includes 
verification   of  past    trends,    and   additional    data   taking   for   better   curve 
definition   at    low  speeds)    is   done   utilizing  only   screen   #9   design.      This 
screen   proved  to  be  most   successful    in   the  prior  work    (References    3,    7, 

10. 

2.      General    Description  of   the   Discharge 

With  no   voltage   applied  across    the   gap,   and  no  gas    flow,    there 
exist   a   small    number  of   random  electrons   present    in    the  gap  which   are 
produced   by  outside   sources    (Reference    12).      As   voltage    is    applied   these 
electrons   under   the    influence  of   the   resultant  electric  field   begin   to 
move   toward  the  anode  producing  a   small    saturation   current  whose  magni- 
tude   is    limited  by   the  number  of  electrons   present   and   their   average 
drift   velocity    in    the   field.      In    a  highly  asymmetric  geometry   such   as 
the  point-to-plane,    the  electric  field  produced    is   very   non-uniform, 
being   highest   at    the  point    surface  and  falling  off   rapidly    in   a   hyper- 
bolic  fashion  with   distance   from   the   point   as    illustrated   by    the  axial 
field  distribution    in    a  hyperboloi dal    point-to-plane  gap    (Reference    17) 
which    is   assumed    in    the   following  analysis: 

EW-^Ep  (I) 

R  =  point  radius 

x=  axial  distance  from  point 

2V 


Ep  =  field  strength  at  point  = 


Rlog(WR) 


V  =  app 1 i  ed  voltage 

d  =  gap  length 
Due  to  this  field  asymmetry,  as  the  voltage  increases  further,  the  field 
strength  at  the  point  soon  reaches  a   value  which  permits  the  free  elec- 
trons moving  toward  the  point  to  gain  enough  energy  to  ionize  by  colli- 
sion some  of  the  gas  atoms  and  molecules  near  the  point  surface,  thus 
additional  electrons  and  positive  ions  are   produced  resulting  in  a  rise 
in  the  current.   The  voltage  at  which  these  first  few  ionization  events 
occur  is  called  the  corona  initiation  voltage  and  is  usually  designated 

V   ,   i.e.,  when  Ep  ~  E   ,   V  s  V   : 
c  c        c 


2V 


Ep  -  Ec  ~  Rlog(WR) 
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E  Rlog(WR) 
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(3) 


From   Reference    13    the   threshold   field     E        can   be   calculated   from   Peek's 
Law   (1929)    modified  by   Gary  et   al    (1972)    for  a   point-to-plane    type 
geometry,    i  .e.  , 


E     =  E  MA 
c  o 
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31    KV/cm      (near   spark   breakdown   of  air) 

coefficient   describing   conductor  surface    (0.6- 

0.308   cm* 

radi  us   of   point 

0.392   P/(273+T)    =    1    for   air   at   STP 
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P   =  pressure    in    torr 
T  =   temperature    i  n      C 


Continued    increases    in   voltage  moves    the  high   field   region    further  and 
further    into   the   gap.      At   the   position  where   the    reduced   field   value 
(E/P)       is   approximately  equal    to   that    required   for    ionization    by  elec- 
tron  collision      (E./P)    ,    free  electrons  will    begin    to  produce   avalanches 
into   the   positive   point.      From   this   critical    distance     x      ,   each 
avalanche  will    produce   an   exponential    increase    in   electrons   and  positive 
ions    left    in    the   aftermath.      The  number  of  electrons   produced    in   the  non- 
uniform field    is   given    by  the  well    known   expression    (Reference    12) 


x2 

/       adx 

xl 
n  =  nQe  (5) 


where 


n      =  original    number  of  electrons 
o  3 

a  =  Town  send 's   first    ionization    coefficient 

x  =  avalanche   travel    distance 

The   coefficient     a      is   a    strong   function   of   the  electric  field  as    illus- 
trated  by   the   usual    semi-analytical    expression    found    in   many    references    ■ 

BP 


a  =  APeEU)  (6) 


where     A     and     B     are   gas   constants,      P      is    pressure.      It    is   necessary, 
therefore,    to    integrate   the  exponent   since    in   the    region    under 
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discussion      E      and    therefore     a     vary  greatly  with     x     as   the  avalanche 
rushes   toward   the  point,    Fig.    3.      An   average  or  effective    ionization 
coefficient   can   be  defined  as    follows: 


a  =  /      a  dx  (7) 

x-x     J 


o 

X 

o 


which  will    produce   the   same   number  of  particles    in    traveling    the   same 
distance,    and   then  we   can   write 


a(x-x   ) 
o 
n   =  n   e 
o 


similar   to    the   simpler  uniform   field   case.      In    the   following   discussion 
the   bar  will    be  dropped   for  convenience,    and     x        assumed  to   be  0.      Due    to 
the  much  greater  electron   drift  velocities,    the  electrons    from  these 
avalanches   are   quickly   drawn    into   the  positive   point    leaving    the   positive 
ions   behind    in    a  maelstrom  of  activity    including  electrostatic   repulsion 
from  the  positive   point,   diffusion    in    all    directions,   neutralization   by 
recombination  with   other  electrons   or   charge  transfer   from  negative    ions, 
etc.      All    of   these  processes   tend   to   dissipate    the   space   charge    concen- 
trations  and   smooth    the  non-uniformity  of   the  electric   field  which    is 
generated  by   the  dense   space  charge  near   the   point,   and   by    the  point 
itself. 

However,    if   the   applied  voltage    is    raised  high   enough    (near 
breakdown),    these  dissipative   processes    lag   behind  production   processes, 


and  the   space   charge   density  and   field   intensity   can   become  of   the  same 
order  of  magnitude   as   the   applied   field.      When   this   happens,    the  portion 
of  the  gap  external    to   the  high   field   region    is   no    longer   able    to  dis- 
tinguish  the   applied  field  nor  the  position  of   the  point,    and   the  dense 
space   charge   actually   assumes    the    role   of   the  point    forming   a  "virtual" 
anode    (Reference    1*0    effectively   diminishing   the  gap  as    seen    from  the 
cathode.      The  effect  of   this   "extension"   of   the  point    is    to    require 
dissipation   of   the  applied  voltage  over  a   shorter  gap,    thus    raising   the 
"drift"    (low  field)    portion   of   the  gap   to   a   uniformly   higher    level    (Fig. 
k) .      Reference    \k    illustrates  how,    in   a  non-uniform  geometry,     the    presence 
of  space   charge   can  modify  an   originally  hyperbolic  field   distribution 
to  produce  a   uniform  field  over  most  of   the  gap    in    concentric   cylinders 
where   the    inner  electrode    radius    is   much    less    than   the  outer. 

Another   important  effect  which  must   be   accounted   for    is    the 
radial    diffusion  of   the  space   charge    in    the    ionization    region  which 
essentially   distorts   the   shape  of   the  original    point,  making   the    resultant 
virtual    anode  appear  extended    in    the   radial    as  well    as    the  axial    dimen- 
sion   (see   Fig.    5) . 

It   will    be   seen    in   the  following   analysis,    that  accounting   for 
the   two  effects    (shortened  gap,    and  modified   radius)    provides   much    in- 
sight   into  various   discharge  phenomena  heretofore   not  well    explained 
and    leads    to   a  very  general    form  of   the  breakdown  equation    for  a  highly 
non-uniform  geometry    in   a    turbulent  flow   field.      The   breakdown  equation 
is    derived  by   accounting   for  modifications   by    the   flow  of   the  quasi- 
static  space  charge  distributions  which   develop    in    the  non-flowing    run- 
ning  discharge   case.      The   trends    revealed  and   results  obtained  are 
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fundamental    and   adaptable   to   different   gas  mixtures,    geometries,    and   flow 
field  parameters. 

3.      Analytical    Development   of    the   Breakdown   Equation   with    Flow 
Interactions 

The  exponential    charged  particle  multiplication   process    in   high 
reduced   fields   has   already   been  mentioned.      Since   the  current    is   produced 
by   the  motion  of   these  particles   between    the  electrodes,    the  current 
density   can   be  written   as   a   summation   of    the   individual    species   contribu- 
tions.     For   the   present    it    is    assumed  that   the  positive    ion   and  electron 
currents   dominate   the  discharge    (Reference   8).      The   current    density  will 
then   be  given   by 

j   =  j .    -  j      =  en.v.    +  en   v  (8) 

■>        J  i        Je  ii  e  e 

where 

j  =  local  current  density 

e  =  electronic  charge  (1.6x10    coulombs) 

v.    =  ion/electron  drift  velocities 
i  ,e 

n.    =  local  charge  densities 
i  ,e 

Considering  only  the  high  field  ionization  region,  if  each  elec- 
tron entering  causes  an  electron  avalanche,  the  average  number  of  elec- 
trons produced  is  equal  to  the  multiplication  factor  plus  one,  while  the 
number  of  positive  ions  produced  would  be  just  the  multiplication  factor. 
The  multiplication  factor  is  usually  so  high,  that  the  number  of  positive 
ions  and  electrons  created  in  this  region  are  identical.   Therefore 

11 


n.  s  n  z   n  (9) 

i    e 


Due  to  the  relatively  low  mobility  of  the  positive  ions,  the 
majority  of  electrons  continue  into  the  anode  and  a  region  of  growing 
positive  space  charge  is  produced  in  the  immediate  vicinity  which  will 
slowly  begin  to  drift  away  due  to  electrostatic  repulsion  from  the  anode. 
Under  these  conditions  the  current  density  can  be  estimated  by  the 
fol lowing 

j  =  en  (vj+ve)  (10) 


Under  conditions  of   low  positive  space   charge,      v        is  much   greater  than 
v.      and   the  equation    reduces    to 


j    ~  env 
J  e 


But  if  the  positive  space  charge  field  becomes  of  the  same  order  of  mag- 
nitude as  the  local  applied  field,  electrons  may  be  retarded  or  even 
trapped  by  the  space  charge  cloud.   In  this  case  the  electron  drift 
velocities  may  be  reduced  to  the  same  order  of  magnitude  as  the  ion 
velocities,  or  to  zero  and  even  negative  if  the  electron  is  trapped  and 
dragged  in  the  opposite  direction  by  high  positive  ion  fields.   In  this 
case,  the  current  density  approaches 

j  ~   env.  (11) 

and  the  particle  density  may  then  vary  between 
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~  _J_ 


n   =  -±-       to     n   s  -J—  (11) 

ev  ev. 

e  i 


It    is   now  proper  to    return    to    the  "anode  extension"   concept    introduced 
above.      As    indicated,   when   the  space   charge   field  near  the   anode 
approaches    the   same  order  of  magnitude  as    the  applied   field,    flux    lines 
from  charges  within  the   remainder  of   the  gap   and   from  the  cathode  will 
no    longer  be  able   to    reach   the  actual    point,   but  will    instead   terminate 
on    charges    in    front  of   the  point,,      This  essentially   shortens    the  gap  and 
to   a   first   approximation,    results    in    dissipation   of   the   applied  potential 
over  a    reduced   distance,    raising    the  mean  electric  field  over   the 
remainder  of   the  gap    (Fig.    5). 

In   order  to   continue  the   development,    an   estimate  of   the    length 
of   the   anode  extension    region    (hereafter  designated  by   a     5    )    is 
required.      A   reasonable  estimate  of   the   beginning  of   the   region  of   dense 
positive   charge  would  be   the  point  at  which    the   applied  electric  field 
approaches    the   breakdown   value    for  the   type  of  gas  mixture    involved.      In 
our  example  of   atmospheric  air,    it    is   assumed  that   the  anode   extension 
in   a  non-flowing   discharge   begins   at   the  value     E(x)    =  25,000  volts/cm 
(E/P    2   35).      This   assures    that   the  majority  of  electrons   approaching 
this   position   avalanche  toward   the  positive   point,   depositing   the   slower 
moving   positive    ions   to   form  the  anode  extension    (see   Fig.    5). 

As  soon  as  a  space  charge  concentration  develops,  however,  there 
are  numerous  forces  and  processes  which  act  on  and  within  it  to  dissipate 
it.      Among   the    important  process    to   be  considered  are: 

(a)      diffusion    in   all    directions, 
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(b)  electrostatic   repulsion   by    the  electrode  of    like  sign,    in 
this   case   the  nearby   point, 

(c)  attraction   by   the  opposite  electrode  which  may   be 
insignificant    in    this   region    in   very  non-uniform  fields, 

(d)  recombination  with  supplementary  electrons  supplied  from 
various  sources,  i.e.,  photoelectrons ,  cathode  emission, 
free  electrons,    detachment   from  negative    ions,   etc. 

In   a   "running"   discharge,    t  .e,  ,   one  which  has    reached  a   pseudo-steady 
mode  of  operation,   a   quasi-equi  1  ibrium  condition    is    ultimately    reached 
in  which    the  charge  production   and  dissipation  mechanisms    locally   balance 
each  other.      The   result  of   this   balance    is   a  quasi-steady   anode  extension 
with   a   finite  and   calculable  extent;   a   distorted   radial    dimension  of   the 
virtual    anode  due   to    radial   diffusion,   and  an    increased  mean   electric 
field    in    the   drift    region  of  the  discharge    (see  Fig.    5).      Reference   3 
alluded   to   these  two   important  dimensional    changes. 

To  calculate   the  equilibrium  boundary  of     6      in    the  axial 
dimension   between   anode  and   cathode,   electrostatic   repulsion,    convection 
(in   the   flowing   case)    and   recombination    must  be   balanced,    diffusion   being 
small    in    the   high   field   region   compared   to   the   axial    component   of   ion 
drift.      The   radial   boundary    is   determined  by   radial    diffusion,    and    re- 
combination,  and    in    the   case  of  crowded  multi-pin  electrodes,    by   symmetry 
consi  derations. 

a.      The  Anode  Extension   -  Axial    Dimension   -   <5 

Since   the    length  of     6      is    determined  primarily   by  a   balance 
between    the    ion   drift  out  of   the  dense  space   charge   region,    and   the  dis- 
sipation  of   that  charge   by    recombination,   an   expression    for   recombination 
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is    required.      The    form  assumed  here   as    in    Reference   3    is    the   three-body 
solution   of   Hinnov  and   Hirschberg,    namely    that    the    recombination    time 


where 


'rec'T^T  (12) 

2pn 


p  =  recombination  coefficient 
n  =  particle  density 
tREr  =  characteristic  recombination  time 

Now,  the  positive  space  charge  created  in  the  vicinity  of  the 
anode  will  be  subject  to  electrostatic  repulsion  until  it  recombines  with 
an  electron.   Therefore,  the  travel  distance  of  an  individual  ion  during 
its  1 i  fet ime  will  be 


xd  =  v.  tREC  (13) 


where      v.    =  k.E   =      ion    drift  velocity.      Inside   this   drift    distance  will 
i  i 

be  much   positive  space   charge,    and  outside    it,   most  of   the   charge  will 
have   been   neutralized  by    recombination.      Therefore,    since 


xd   =  5  (U) 


comb  in  ing  Equations    12,    13,   and    \k 


5  =   v,(-!-jO  (15) 

2pn 


Substituting    for     n      from  Equation    11 


v.    ev     2  v.    ev.    2 

6  =  xU-A        to     M-J4  (16) 

2p     J  2p     j 


So  far,  flow  has  not  been  introduced  so  the  above  equation  is 
general.   Taking  flow  into  account,  the  following  notation  will  be  used 


6   =  anode  extension  in  flowing  gas 

5   =  anode  extension  with  no  flow 
o 

5 

5*  =   flow  value  normalized  by  no  flow  value   (■?— ) 

o 

Typical    flow  velocities    in   practice   range   up    to   about      1-2x10      cm/sec. 
Typical    ion   drift   velocities    in   atmospheric  air  at   the   outer    limits   of 
the   anode  extension      (E/P  =  25,000/760   ~  33)    also   run    in   the  neighborhood 
of  5-8x10      cm/sec. 


Defining    the   slip   parameter     S   =  u/v.    =  r— p-     ,    significant   flow 

i 
effects   can   be  expected  when    the   slip   parameter    is   of  order  one,   especial- 
ly   in   a   dense   collision   dominated   discharge  where    the  mean   neutral    gas 
motion    is  easily   transfered    to   the    ion    population.      As   can   be    seen    in    the 
speci  f  i  c  examp  le  , 


S   =  JL_::J£!2£-  2x10^    il3   .|4J,  (17) 

k.E        8x10  5x10 


and  substantial    space   charge   redistribution   should  be  anticipated. 

Going   back   to   Equation    16,   and    rewriting    to    include   flow  effects, 
we  have 

a.      No  electron/ion    coupling: 
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and 


(v.+u)      ,   (  k.E(l+S)v  2 

6u  =  "V-e    %)  "      '     .    2  %>  (18) 

J  J 

Ju  J  u 


v.v  (k.E)v 


5*.  JL-U4iL  (18b) 

o  j* 


b.      Heavy  electron/ion    coupling    (high   space   charge   fields) 


(v,^)3  (k.E)3(1^S)3(^_)  , 

u  .    2      V2D;  ^1  V2p'  U3' 

J  J 


J  J 


6^  =  iH|L  (19b) 

j* 


where  J  ,s  =  ■: — 

Jo 


From  these  considerations  it  can  be  seen  that  the  power  of  the 
velocity  term   (1+S)   will  vary  with  space  charge  density  (or  current) 
and  we  will  write  a  generalization  of  1 8b  and  19b  as  follows: 


5*  =  i!±sT  (20) 

j* 
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and 


where  m  is  a  variable  exponent. 

Since  continuity  considerations  require  that  the  total  current 
passing  any  equipotent ial  surface  within  the  discharge  be  constant  the 
total  current  can  be  represented  by  the  electron  current  entering  the 
anode  where  ion  contributions  are  zero,  so 


Xc  -   .  af(<5)6  f-,\ 

I  =  I  e        I  e  (21 ) 

o         o 


where  j    is  the  electron  current  flowing  into  the  ionization  region 
(defined  by    x    which  is  proportional  to  5  )  ,  the  exponential  is 
the  multiplication  factor,  and  f   is  an  as  yet  unspecified  function  of 
5  ,  therefore  - 


a  f<5 

with   flow:  j      =  j      e  (22a) 

Ju       Jou 

aof5o 
no   flow:  j      =   j      e  (22b) 

Jo       Joo 


fa  5   (a*5*-l) 
then  j*  =  jQ*e     °  °  (23) 


NOTE:   The  *   value  of  a  quantity  is  used  consistently  to  indicate  the 

value  of  that  quantity  in  a  flow  field  normalized  by  its  no 

flow  value  (except  for  pressure  -  P*  to  be  introduced  later  is 

equal  to  P/P  =  local  value/value  at  STP  =  760  TORR) . 
o 
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The   function      f(5)      has   been    introduced  to    relate   the   critical 

ionization   distance    (that   point   at     E(x)      =     EAppL| EQ  +   EspACE    CHARQE     S 

E.)      to    the   anode  extension     <5    .       It    is   apparent    that     6      is   proportional 

to  x  i.e.,   as     <S      increases  either  due   to   convection   or    increasing 

c 

applied  voltage,        x  also    increases   due   to    increased  applied  or  dis- 

placed  space   charge   field,   and   conversely,   as    the   anode  extension   shrinks, 

x         will    decrease,,      However,    due    to   the  non-linearity  of   the  fields 
c  '  ' 

involved,    a  variation    in     5      will    not   produce   a    linear   change    in      x 
To  account    for   the  non-linearity,      f(<5)       is   assumed   to  obey   a   simple 
power   law   form 

f   =   CS*n  (2k) 

Squaring     j*     and   substituting    into    Equation   20,    yields    the   following: 


5* U±52? (25) 

Jo"Zexp[(2CS"nao«o(cx*4*-l)] 


Rearranging   and   taking    the  natural    log  of  both   sides, 


£n<5*  +  2Ca  5  <5*n(a*6*-l)   =  mfcn(1+S)    -  In}    ~'*1 
o  o  o 


and  solving  for  5'    and  collecting  6*  terms  on  the  left  side, 


S""   -  S~^r  *   <2*rV-)i*  »"«.  -  (^_)>(,+S)  -  *nJ0*2     (26) 
0  0  o  o 


15 


This  is  the  basic  anode  extension  equation  in  a  flow  field,  where  again 
the  *  values  are  the  flow  values  compared  to  the  no  flow  values,  i.e.: 


6*  =  5  /6    =  normalized  anode  extension 
u  o 

a*  =  a  /a    =  normalized  ionization  coefficient 
u  o 


S  =  u/v.    =  the  ion  slip  parameter 

i  *  =  J   /J   -  normalized  saturation  current 
Jo    Jou  Joo 

C  =  CONSTANT 

The  exponent  n  of  6*  turns  out  to  be  a  free  parameter  of 
the  problem  whose  exact  nature  is  not  yet  well  established,  but  is  under 
investigation.   The  parameter  m  is  the  power  of  the  velocity  term 
(1+S)   which  analytically  was  determined  to  range  from  approximately   1 
to  3   ,  although  experimental  data  indicate  that  occasionally  a  higher 
value  is  required.   The  reasons  for  this  disparity  are   also  under  active 
investigation. 

Inspiteof  the  limitations  cited  above,  this  equation,  even  in 
its  preliminary  form  (Equation  26)  describes  all  the  previously  observed 
breakdown  voltage  characteristics  on  the  particular  apparatus  used  in 
the  initial  tests  for  various  velocities,  turbulence  levels,  and  flow  ori 
entations  tested,  even  including  some  seemingly  anomalous  points  in 
earlier  tests  which  radically  departed  from  the  norm  as  reported  in 
References  3  and  7.  The  equation  also  predicts  some  salient  features  in 
the  low  velocity  ranges  which  had  not  previously  been  observed  due  to 
strong  emphasis  on  high  speed  testing.   These  low  velocity  features  have 
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been    confirmed    in    recent    tests   and  will    be  described    in   detail    in   part     6 
of   this    section  of   the   report. 

To   check   the  mathematical    correctness  of    the   form  of   this   equa- 
tion,  when   there    is   no   flow,      S  =   0    ,    and   by   definition,    the    (*)    terms   = 
1    .      Therefore,    the   three     Jin      terms   go   to  zero,        a*  =    1      and  we  are 
left  with   the    identity      0=0     as    it    should. 

The   peculiar  geometric  factors   of  a   problem  are   all    contained    in 

the     6      terms.      Gas   and  electric  field  parameters  enter   through    the     a 

and     S      terms;   and  turbulence  effects  enter   through     a*     as  will    be 

explained    in    the   following   sections.      Secondary  electron   production 

processes   are    included    in      I    *  . 

o 

To  obtain   a   solution   to   Equation   26,    requires   some   auxilliary 

equations  which   are   developed  below. 

h.      Auxilliary  Equations 

a.      The  no-flow  anode  extension      (<5    ) : 

o 

The  basic  assumption  of  the  anode  extension  model  is  that 
the  gap  breakdown  potential  is  a  linear  function  of  the  effective  gap 
length,  i.e.,  the  actual  gap  length  less  the  anode  extension  (which  in 
effect  is  a  portion  of  the  gap  already  in  a  state  of  local  breakdown). 
Therefore: 


Vb  =  d"  z   d  -  o  =  E(d-5  )  +  V  (27) 

o    ■        o        o     c 


Vb  =  no  flow  breakdown  voltage 


d'  =  effective  gap 
d  =  actual  gap  length 
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<S  =  no  flow  anode  extension 
o 


E  ~  mean  electric  field  in  effective  gap 
V  5  corona  initiation  voltage 


For  a  specific  geometry,   Vb   can  be  calculated  from  the  condition  that 

E(x)  =  E.   at  x  =  C'S      ,  resulting  in 
i  o  '        3 


cf    \        E.  d        E  d  ,        2Vb 

E(X)     ~         I     „  R  p  R  1  O  /,n\ 

~~P        ~  V  "  TfHxT  P  (R+C6    )    7  Rlog(WR)  K      ' 


and 


VbQ  =  IlP*(R+C-5o)log(4d/R)  (29) 


for   the  hyperboloi dal    point-to-plane. 

An   analytical    form  can    also  be   determined   for     V        as    in 
Section   2   by   assuming    that      V        is   a    linear   function  of   the  point   field, 
i.e.,   when    the    region    in    the  vicinity   of   the   point    reaches   the   reduced 
field   strength    for   ionization.      Again    for   the  hyperboloi dal    point,    this 
results    in 


E  2V  E. 

■P2  =    PRWWRi   "    Cl  F*  (30) 


and 


C   E. 

VC  ~  4~Lp,VR1og(ltd/R)  (31) 


but  Vb      -   V     =   E(d-6    ) 

o  c  o 
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and    inserting   Equations   29   and  31    into  27  yields 


<5     z 
o 


.  E. 
C      i 


P*(j~!-)log(J»d/R)    +    1 


(32) 


b.       Ionization    Coefficient   Equations    (a  and  a*) : 

(1)      No   Flow  -   use   Equation  6  with      0     subscript 


a     =  APe 
o 


Bp 
1 


(33) 


(2)      With   Flow/Turbulence   Effects     a      : 

u 

The   primary   coupling  of     a     with    the    flow   is    result   of 

pressure/density   fluctuations    in   a   highly   turbulent   flow.      Therefore,    by 

substituting      P  ±  AP      into   the   above  equation   and    taking   a   simple   average 

a(P+AP)   +  a(P-AP)j  /2    ,    and  eliminating   second  order   terms    results    in    - 


a.     ~  tt 
u       2 


BAP  BAP 

E     a        E 

i  +  e 


(34) 


or 


a 


BAP 
E 


BAP 


+  e 


(35) 


The   pressure    fluctuations   due   to   turbulence   can   be    related   to   the  mean 
flow  velocity    through    the    local    pressure-density    relation    in   an    incom- 
pressible  flow,    i.e.,    the   total    pressure  equals    the   sum  of   the    static 
and   dynamic  pressures. 
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1        2 

p  ■  p0  +  r  pu 


(36) 


Assuming    local    test    section    total    pressure    remains    constant    in    the 
presence  of   turbulence,    and    that     vu'    ,    the   fluctuating   component  of 
velocity,    is   proportional    to     U    ,    the  mean    flow  velocity,    then 


P   -   P     =  AP   ^   £pu'2 
o 


(37) 


The  quantity      (vu"   /U)      is    called  the   turbulence    intensity    (Reference 
18).      At   constant   temperature,      p   ~   P    ;    therefore   for   constant  T    , 


AP   *   ?pU2 


(38) 


but   also 


S  ■   U/k.E  +  U  =  Sk.E 
i  i 


u/k.E     f     u     =     U  -   Sk.E 
i  i 


k.    =   k     «SL 

i  o  P 


and 


AP  =  EP(Sk.E) 


(39) 


Substituting    into  Equation   35   yields    finally   - 


a*  =  tt 


S2E 

B?poko  hr  + 

e  +e 


2  2 
B£P  k      S   E 
o  o 


p* 


(AO) 


which  when      u   =  U   =  S   =   0    ,    reduces    to     a-  =    1      as    it   should. 
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c.   Breakdown  Voltage  Equation  with  Flow  and  Turbulence 

The  ultimate  objective  of  the  above  derivations  is  to  arrive 
at  an  expression  for  computing  the  fundamental  gas  dynamic  effects  on  the 
breakdown  voltage  of  an  electric  discharge.   Approximating  the  breakdown 
voltage  ratio  as 


Vb    d-6 

vb* =  vet  ~  cRT  from  Equation  2?  »  (^0 

o  o 


then    in    final    form, 


6 
1   -   (/)6* 

Vb*  =  b (42) 

o 

(i  -/) 


which   also   reduces    to     Vb-    =  1      when      U  =  0  ■+  6*  =    1    . 

It  should  also  be  noted   that  when     5*  -*■  0,    Vb*)  ■»■  «—  ,      the 

max  o 

same    relation   which  was   postulated    in    Reference   3   based  on    the   anode 
extension    concept. 

d.      Summary  of  Equations 

For  convenience  of   reference    in   the   following   discussion, 
the   basic  breakdown   potential    and  auxilliary  equations   are    repeated 
below   : 
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(1)   Normalized  Anode  Extension  <5*  =  ■*—  : 


i*"  -  ^^-+  (fcsV'^4  -  <2toV)>(,+s)  -  inv2        <2S> 


o  o 


o  o 


(2)      No   Flow  Anode   Extension      6 


5     =  tr 

o  „  t . 


P*(|.J.)log(WR)    + 


(32) 


(3)   No  Flow  Ionization  Coefficient  a 


a     =  PPe 
o 


BP 
E 


(33) 


a 

(k)      Normalized    Ionization    Coefficient     a*=  — 

a 
o 


a*  =  — 

a    2 


2  2  2  2 

B£P  k  S  E    -B£P  k  S  E 
o  o  o  o 


(40) 


Vb 


(5)   Normalized  Breakdown  Voltage  Vb*  =  -rrr— 


Vb 


o  ** 


Vb*  = 


6 


(M) 
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5.   Qualitative  Discussion  of  the  Equations 

The  key  to  a  clear  understanding  of  the  equations  lies  in  a 
simple  analysis  of  the  normalized  anode  extension  6*  .   The  basic  form 
of  this  equation  is 


«*  -  f(|r)  C3) 


In    fact    it  was   derived   from  Equation   20. 


2  2 

But     S  ~    u     and     j*     **   n*       which  gives 


m 


which    is   essentially   the   relationship   derived    in    Reference   3.      Thus,    the 
anode  extension    is    a   product  of  two   factors,    the   pure   convective   velocity 
which,    in   normal    flow,    tends    to    lengthen      8*    ,    and  the    inverse  of   the 
charged  particle   density,   which    tends    to   shorten     6*     due   to   faster 
recombination   at  higher     n*     and   consequently    results    in    less   particle 
dri  ft. 

In    the    reverse    flow  mode,    the     <5*     equation    can   be   approximated 
by: 

6*  =    (1+S)m  (45) 
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where  u  and  therefore  S  are  negative,  which  decreases  5*  in  a  rapid 
fashion  with  increasing  negative  velocity.   It  is  interesting  to  note  that 
in  the  reverse  flow  approximation,  a*  cancels  out  and  Vb*  is  not 
influenced  by  turbulence. 

The  breakdown  voltage  behavior  is  opposite  in  sense  to  5*  ,  so 
when  8*      increases,  Vb*  decreases  and  vice-versa,,   Note:   frequent 
reference  to  Figure  6  wi  1  1  be  helpful  in  interpreting  the  following.   As 
mentioned  in  the  reverse  flow  mode,  8*      is  quickly  reduced  and  therefore, 
Vb*  immediately  rises  (Part  2  of  Figure  6)  and  quickly  approaches  an 
asymptotic  limit  (Part  1  of  Figure  6)   as  8*  •*   0  corresponding  to  the 
requirement  to  breakdown  the  entire  physical  gap. 

In  the  normal  flow  mode  (flow  from  anode  to  cathode),  the  initial 
low  gas  velocity  tends  to  lengthen  8*    ,  which  causes  a  decrease  in  break- 
down voltage  (Part  k   of  Figure  6),  but  due  to  the  exponential  dependence 
of  n*  ,  j*  ,  and  a*  on  U   (through  S)  ,  with  effective  turbulence  the 
a*  terms  in  the  denominator  of  Equation  26  soon  prevail,  causing  6*  to 
decrease  and  Vb*  to  increase  at  higher  velocities  and  turbulence  levels 
(Part  5  of  Figure  6).   The  velocity  at  which  this  transition  occurs  is  a 
direct  function  of  the  effectiveness  of  the  turbulence  introduced.   If  the 
turbulence  is  inhomogeneous  or  ineffective,  the  Vb*  curve  may  not  reach 
a  minimum  in  which  case,  further  increases  in  velocity  may  actually  reduce 
power  input  to  the  discharge  (Part  7  of  Figure  6,  dashed  curve  extension). 

It  is  interesting  to  note  also  that  in  the  normal  flow  mode  with 
effective  turbulence,  after  initially  increasing,  5-  can  again  approach 
0   as  its  absolute  limit  (Part  6  of  Figure  6)  due  to  the  high  recombination 
rate,  thereby  confining  the  normal  flow  Vb*  rise  to  the  same  asymtotic 
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limit  as  applies  in  the  reverse  flow  case,  and  again  corresponding  to  the 
requirement  to  breakdown  nearly  the  entire  physical  gap. 

It  should  also  be  noted  that  this  asymptotic  behavior,  which  has 
been  observed  in  many  cases,  (References  2  and  3),  poses  a  fundamental 
limitation  on  the  stabilization  of  the  discharge  potential  by  gas  dynamic 
means.  The  final  breakdown  may  possibly  be  the  result  of  either  of  two 
separate  sequences: 

(a)  the  relatively  uniform  external  "drift"  field  may  reach 
such  a  high  value  that  a  random  electron  avalanche  originating  in  the 
vicinity  of  the  cathode  can  take  enough  energy  from  the  field  to  com- 
pletely cross  the  interelectrode  gap,  resulting  in  an  internal  break- 
down of  the  gas,  or, 

(b)  if  the  cathode  material  has  a  low  work  function  it  may 
begin  emitting  numerous  secondary  electrons  prior  to  the  drift  field 
reaching  the  gas  breakdown  potential  in  which  case  the  breakdown  proba- 
bility will  rise  and  a  cathode  limited  breakdown  may  occur. 

As  turbulence  intensities  are  increased  for  a  given  gap  length, 
the  reverse  flow  region  is  virtually  unaffected,  except  to  hasten  the 
rise  of  Vb*  to  the  asymptote,  while  in  the  normal  flow  direction, 
increasing  turbulence  results  in  higher  a*  values  with  an  eventual 
upswing  in  Vb*  due  to  the  resultant  6*      reduction.   Figure  7  illus- 
trates the  progressive  modification  of  the  Vb*  curves  with  increasing 
turbulence.   The  asymptotic  behavior  observed  at  high  normal  velocities 
and  turbulence  levels  confirms  the  apparent  rise  in  the  Vb*  curve  and 
explains  logically  the  apparently  anomalous  "short  gap,  high  current" 
data  points  reported  in  References  3  and  7. 
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The   analytic  solutions    illustrated    in   Figure   7  also  provided   an 
additional    area   for   test  and  verification  of   the  breakdown  equations.      It 
can   be   seen   from   the   Figure   that  even    in   the  "short  gap,    high    current" 
cases,    there  should   be  at    least   a   small    region  of  decreasing      Vb*     at  very 
low  normal    flow  velocities  due  to   the  non-turbulent   convective  effect. 
Recent    low  speed   tests   have   confirmed   this   characteristic  as    the   data    in 
Figure   8   show. 

Another  qualitative  assessment  of   the  equations   can   be  made   by 
observing    the  behavior  of   the  curves   as   a   function  of   the  gap    length. 
Figure  9    illustrates    this   behavior   for  a    1    and  5   centimeter  gap  which 
encompasses    the   ranges    tested.      The   trends  of  the   curve   shift    in   both 
normal    and    reverse    flow  with   changing  gap    length  were   as    indicated   by    the 
theoretical    curves. 

Section  6   provides   a  more  quantitative   comparison   of  observed 
data  and   the  appropriate   theoretical    curves.      However,    before  moving    into 
that  subject,    a   very    interesting  and    rather  striking  qualitative   comparison 
can   be  made  of   the  form  of   Equation   26   and  Meek's   Breakdown   Criterion 
(as   presented    in    Reference    12)    for   the   uniform   field    in    a   static    (non- 
flowing)    discharge»      As  will    be   seen,   despite   the   very  different    conditions 
in    the   two  discharge   situations,    and   the  very  different    approaches   used   in 
deriving   the   breakdown  criteria,   the    resulting   breakdown  equations   are 
remarkably  similar. 

Since  a  description  of   the  non-uniform   field  anode  extension   has 
already   been   developed,    a   short  description   of   the  Meek-breakdown   scenario 
in   the  uni  form   f iel d  wi 1 1    help  put    the   remainder  of   thi  s    section    in 
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perspective.      The  Meek   Criterion    is   based  on    the   relative   space    charge 
field   strength   at    the  head  of  an   electron   avalanche   compared   to   the  applied 
electric   field.       If   the  avalanche  travels   a   certain   critical    distance   under 
the    influence  of  the  applied   field,    the  space   charge   field  at   the  surface 
of  the  space   charge    region  will    become   the  same  order  of  magnitude   as    the 
applied   field.      When   the   field   ratio     K  =  Er/E  "^  1-1,    corresponding   to    the 
avalanche   traveling    the   critical    distance     x      ,    the   space   charge   field   can 
no    longer  be    ignored.      The   Meek   Criterion    from  Reference    12    is  written   as 
fol lows : 

ax  +  in-  =    ]k.kG  +  Jin^-  -0.5£npx  +  in   x  (46) 

p  pd 

The   method  of   solution   and   some    results    are   discussed    in    detail 
in   the    reference.      What  we  will    concentrate  on   are    some   straight   forward 
substitutions    into  the   above  equation    for   comparison   with    the  anode 
extension   model.       If    the   gap    length    is   adjusted   such    that   an   avalanche 
becomes    critical    just   as    it    reaches    the   anode,    then    the  electric   field 
lines  will    be   distorted    in    the   direction  of   the   space   charge,    as    in 
Figure    10a   and   the   concentration   of   charge  will    produce   a  non-uniform 
field  and,    in   fact,    just   as    the  avalanche   touches    the   anode,    it  acts   as   a 
small    area    (like  a   point)    of  anode  extension.      The  axial    length  of   this 
extension    is      2r    ,    but      r      (the  avalanche    radius)    is   proportional    to   the 
distance    it  has    traveled,    so   from     2r  *"   5   ~   x    , 

x  =   cS  (47) 
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and    it   can   already   be   seen    that  the  situation    looks   very  much    like  the 
point-to-plane  geometry.      After  making   this   substitution,    and   then 
rearranging   and  substituting  non-dimensional i zed   flow  variables    for  all 
the  quantitites,    the   result    is 

5*  -    (      l—jJUnfi*  =  -(_*     )i-£ns  +    (—V-)-^  "    (-^)-Lm  °       %      (W) 

ca  6     a»  ca  5     a-  ca  6     a-         a  6     a*  u 

oo  oo  oo  oo 

-    compare  with   Equation   26   for  the   non-uniform   flow   field   - 


o  o  o  o 


The    transformed  Meek   Criterion    looks   very  much    like  equation   26  with    some 
sign   changes,    but   the   basic   form    is  nearly    identical    to   the  non-uniform 
flow  field  equation.      The   sign   changes   can   be   accounted   for   by   the  obser- 
vation  that    in    the   non-uniform   field  equation: 


6o  ~    logCtd/R)  (49) 


and  as  the  point  radius  (sphere  or  wire  to  plane  equations  are  similar) 
R  ■*■  °°  ,  the  log  term  becomes  negative  and  the  equations  merge  to  nearly 
the  same  form. 

After  seeing  how  similar  the  actual  physical  situations  are  at 
breakdown  though,  the  similarity  should  not  come  as  a  surprise.   In  fact, 
what  the  similarity  further  suggests  is  that  perhaps  there  really  is  no 
absolute  uniform  field  breakdown.   In  a  uniform  gap  at  a  potential  just 
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below  breakdown,    statistically,   numerous   small   electron    avalanches  may 
originate  and   die  out    in   a   short   distance   as    depicted   in    Figure    10a.      Many 
of   these  mini -avalanches  must  meet    the  anode   to  form  minute   deformations 
of  the   field   there,   any  one  of  which,    if    large  enough,    can    ultimately 
trigger  total   gap   breakdown   similar   to   the  non-uniform  point-plane. 

This   brings   us   to   an    interesting  obeservation    concerning   a   highly 
non-uniform   field   situation  where   the  accumulation  of  space   charge   around 
a  positive  point   tends    to  distort    its   apparent   shape  and    radius,    the   result 
of  which    is   to   tend   to   reduce   the  non-uniformity  and  go   toward   a  more    uni- 
form field   configuration.      From  these  two  examples,    it  would   appear    that 
the  "natural"   state   at   breakdown    in   a   discharge    is   a  moderate   state   of 
non-uniformity   rather   than   at  either  extreme. 

6.      Quantitative   Comparison  of   Data   and  Theory: 

Figures    11    through    17  present   some    typical    data    recorded  on 
various    runs  with   the   point   and  wire-to-plane  electrodes.      Conditions 
ranged   from  high    reverse   flow  velocities      (-120  m/s)    to  high   normal    flow 
(80-100  m/s  with   turbulence).      No   turbulence   could  be    introduced    in    the 
reverse   flow  mode,    since   as    illustrated    in    Fig,    1,    the   cathode    is    constructed 
of  aluminum  honeycomb   flow  straightener  which   assures    smooth    laminar   reverse 
flow   throughout    the  velocity   range   tested. 

The  data   presented  were   collected  over  a   six  week   period    in 
varying  weather   conditions   and  are  not   corrected   for   relative   humidity   and 
ambient    temperature  effects,    both  of  which   affect   the   properties    being 
measured.      Correction   for   these  effects  might    reduce   the   scatter  evident 
in   various   data   sets,   especially   in    the   composite   summary   plots.      An 
additional    source  of  statistical    scatter  due   to   unequal    current    flow   to 
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various   pins   or  wires   may  also  be   significant,   especially    in   a   discharge 
such   as  ours  which    is   not   ballasted  to  produce   uniform  current   flow   from 
each  element.      As  noted    in    Ref.    3,    however,   effective   turbulence   does 
cause  more   uniform  activation   of  all    pins   and   therefore    reduces   the 
adverse      effects   due   to  no  ballast. 

Figure    11    displays    corona    initiation    and  breakdown   voltages   with 
no   flow  versus   gap    length   for  both   geometries   tested.      Theoretical    curves 
are  also  displayed  which  were  computed   from  the   appropriate  equations 
developed    in   Section  A      .      These   curves   display  an    intersection   of   the 
Vb     and     Vc     characteristics   at   the   shorter  gap    lengths   and  are   very 
similar    in   appearance   to  those  presented    in    Reference   8   for  single 
point-plane  geometries.      The    intersections   define  what    is   termed   the 
"existence"    region    in    Reference   8   for   the   corona   discharge*,      To   the 
left  of   this   point,    a  corona  discharge  will   not    form,    and   the  gap  will 
breakdown    (arc)    upon    reaching   the  value  of     Vc   .      To    the   right   of   that 
point   a   corona  will    form  which  enhances    the  discharge   current.      A  common 
empirical   expression  which   has   been   found   to  quite   closely    fit    current 
data   for  various  non-uniform  geometries   can   be  written   as    follows: 

I    =  KCV(V-Vc) 

The    initial    constants   contain  geometric  and  gas   parameters.      This   equation 
highlights   the   dependence  of    I    on   V  and  on    the  difference   between    the 
applied  and  corona    initiation  voltages.      The  existence    regions    for   the 
present   configurations   are  seen    to  be   for  - 

Points-Plane        d  >  0.2      cm 

Wi  res-Plane  d   >  2.0   cm. 
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Figure  12  presents  a  comparison  of  Vb*)     from  Equation  hi 

max 

when  5*  -*■  0  .   Due  to  power  supply  limitations,  insufficient  points  were 
accessible  in  the  wire/plane  configuration  because  of  higher  breakdown 
voltage  requirements  to  adequately  define  the  wire/plane  curve  so  only  the 
point/plane  data  are  presented.   This  curve  is  interesting  in  that  it 
shows  at  longer  gaps  flow  cannot  affect  breakdown  voltage  as  radically  as 
it  does  in  a  short  gap.   The  general  trends  of  the  Vb*  vs.   velocity 
data  taken  under  similar  conditions  of  humidity  and  temperature  tend  to 
confirm  this.  However  some  contradictions  are   evident  and  this  cannot 
be  considered  a  conclusive  result. 

Figures  13,  1^,  and  15  are  comparisons  of  data  and  theory  for 
three  gap  lengths  of  Vb*  vs.   velocity  .  The  values  of  the  parameters 
C,  m,  n,  and  turbulence  index  used  in  the  calculated  curves  are  given.   It 
has  been  observed  that  C  =  1.3  ,  m  =  3  ,  and  n=  8  generally  result  in 
a  good  first  estimate  for  the  basic  curve,  and  then  the  turbulence  index 
(input  of  turbulence  intensity)  is  varied  to  match  the  data.   In  using 
this  technique,  it  has  been  observed  that  long  gaps  usually  require  a 
lower  turbulence  index  even  though  the  actual  turbulence  being  generated 
at  a  given  velocity  is  the  same  as  in  shorter  gaps.   One  possible  explana- 
tion for  this  is  that,  as  pointed  out  by  Reference  2  after  the  turbulencee 
being  generated  is  fully  developed,  it  begins  to  rapidly  decay  with 
time  and  distance,  so  in  a  long  gap,  although  turbulence  levels  entering 
the  discharge  region  are  the  same  as  in  short  gaps,  the  turbulence 
level  leaving  the  discharge  will  have  decayed  more,  resulting  in  a  lower 
average  turbulence  intensity  in  the  longer  gap.  This  would  indicate  that 
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a  weighting  of  the  turbulence  index  with  gap  length  may  be  required. 
Probe  measurements  with  screen  #9  installed  in  the  test  tunnel  show  that 
at  the  highest  attainable  velocities,  the  turbulence  intensity  peaks  just 
at  the  beginning  of  the  discharge  region  (at  the  pin  tips),  which  would 
support  the  above  argument,  and  which  also  may  be  the  key  to  the  success 
of  the  screen  #9  design. 

Figures  16  a,  b  and  17  a  and  b  are  summary  plots  of  Vb  and   lb 
data  vs.   velocity  for  the  points  and  wires.  While  many  gap  lengths  and 
turbulence  levels  are  included  under  varying  ambient  conditions,  the  plots 
serve  to  point  out  the  main  features  and  general  trends  of  these  parameters 
with  flow  and  turbulence.   The  lines  on  these  curves  are   not  calculated, 
but  are  drawn  in  to  highlight  the  basic  features  of  the  various  data 
collected  and  to  show  the  similarity  with  the  theoretical  curves  dis- 
played elsewhere. 

The  next  section  qualitatively  examines  the  wire/plane  in  more 
de  ta  i  1 . 
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B.      Wi  re-to-P  lane  Geometry 

1 .      Qualitative   Discussion  of   Differences    from  Point-Plane 

One  of  the  significant  changes    in   going   from   the   point-plane 
geometry    to   the  wire-plane    is   a  major  modification    in    the  shape   of    the 
electric   field  strength  distribution   across   a   gap  of  given    dimension   and 
a   fixed  applied  voltage.     While   the  assumed  distribution    (proportional    to 
x     )    is   still    "hyperbolic"   as    illustrated   in    Figure   3,    at   the  electrode 
surfaces    the   point    field   is  much   higher  and   drops  off  much   more   quickly 
than    a  wire    field   of    the   same    radius   at    the   same   potential.      Since    in 
both   cases   the   cathode    is   grounded,    the  wire  electric   field    in    the   drift 
region   must   be   higher   than    that   of   the  point-plane.      The  wire    radii    used 
in    the   current    research    started  at  0.08   cm,    but    it  was    found   that    to 
operate  within    power  supply   capabilities,    the  wire   size   had    to   be    reduced. 
As   a   result   the  wires  were    incrementally   filed   down   to  a  wedge   shape  of 
approximately    30      angle   and  maximum  thickness   of   about   0.05    cm.      This 
permitted   some   current  enhancement  within    available   power  supply    limits. 

Several    things   are  evident    in    comparing  point  and  wire   distribu- 
tions   (see   Fig.    3   ).      First,    for   a  given   gas,    the   ionization    field   strength 

E. determines    the    length  of   the  anode  extension     6    ,    so    i f     E.         is    lower 
i  i2 

than    the  point  of   crossing    (point   P)    of   the   curves   at  a  given    voltage, 

6       will    be  greater   than     6      ,    but    if     E.         is   higher   than    this    point, 
w  a  p    '  i  ,  3 

then     <5       will    be  qreater  than     <5      .      Since  as  will    be   discussed    later, 
p  w 

the   current    is   a    function  of     <5    ,    this   behavior  can   have    important    con- 
sequences. 

Another   important   characteristic  difference  of   the   point    and  wire 
distributions    lies    in    the   transition    region    from  the  high    field   to   the  drift 
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field.      The   transition    is  well    defined   for  a  sharp  point,    but   becomes 
very  slow  and    indefinite   in  wires   and  points   of   large   radii.      This 
results    in   a    less  well    defined  anode  extension    in    the  flatter  distributions, 
which  should  affect   the   flow   interactions    in    the   loosely   bound  outer  por- 
tions of   the  space  charge      (<5)      region.      Changes    in    the  behavior  of   the     6 
region  with   flow  should  show  up   as   changes    in      Vb*     with  velocity.      Some 
major  changes    in    the  wire-plane  and  point-plane   Vb*  vs.    velocity   curves 
have  been  well   documented   in   the   current    research.      Figure    18   illustrates 
the   basic  nature  of   these  changes. 

It  has   been   consistently  noted   that   the  wire   curves,   while   still 
exhibiting   the    low  velocity  depression   of     Vb*   ,    almost   always    take   an 
early    (i.e.,    low  velocity)    upturn    toward   the   asymptote  whether   the   flow   is 
turbulent  or  not    (Fig.    8).      The   convective   depression,   caused  by   an    initial 
lengthening  of     5    ,    is   very   small,   and   it    is   thought   that   the    rapid    rise   in 
Vb"      is   due   to  a    rather    large  movement   and  dissipation   of   loose   space 
charge  at   the  outer   limits  of   the  anode  extension.      In    this    region,    due   to 
the   flatness   of  the  E-field,    the   ion   slip  parameter      (S  =   u/k.F.)      can   attain 
a   significant  value  over  much  of     5      at  once,  with   the    result    that    the  outer 
portion  of   the   anode  extension   actually  blows  off  or  becomes    too   tenuous   at 
low  velocity   to  act  as  an  extension,  which   then   causes   a   large   and   rapid   rise 
of     Vb*   .      The   asymptotic   limits  of     Vb*      in   both   the  normal    and    reverse 
flow  directions  are   interpreted  again  as   the   requirement    to  break  down   the 
entire  physical   gap,   as    in   the  point-to-plane   development. 

A  third  characteristic  of   the  wire  geometry    is    the    inherently 
higher  values  of  electric  field   in    the  drift  portion   of   the  gap.      This 
increase    in   the  drift  field    increases    the  probability   that  high    ionization 
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and   streamer   activity    in    the   region  near  the  wire  will    be  able    to  develop 
and  sustain   self-propagation   through   the  drift   field   region    to   the   cathode 
by    the   streamer  mechanism  discussed    in  most   texts  on    the   subject   (see   Refer- 
ences  8,    12,    15»    16)    resulting    in   breakdown.      On    the  other  hand,    however,    in 
a   discharge  designed   for  optimum  energy   transfer   to  a    laser  gas    (pumping), 
the  design  of  the  electrodes  will    to  some  extent   be  dictated   by    this   pumping 
requirement.      In   practice,    a   trade-off   is    required  between  electrode   design 
optimized  for  total    power   input    to   the  gas   and   that   design   optimized   for 
best   E/N    in    the   drift    region    for  most  efficient   pumping.      In   an  electric 
discharge    laser  cavity,    for   instance,    an    input   power  of    1000  watts/cm 
means   very    little    if  only   two  percent    is   going   into   the    lasing   states    due 
to  an    inefficient  E/N    in   the  pumping    region. 

One  of   the   attractive   features   of   the  wire-to-plane  geometry    is 
the   increased   ionization   volume  available   continuously   along    the  wire.      The 
point-to-plane    ionization   volume,   on    the  other   hand,    is   confined   to  what 
would  amount   to  discrete  points   along    the  wire. 

The   point-plane  produces   discrete   3~dimensional    axial ly   symmetric 
multiplication   volumes,    at  each   point,   while   the  wire-plane   system    is   a  2- 
dimensional    continuous   volume   spanning   the    length  of   the  wires.    Most   elec- 
tron multiplication  occurs  within   these  volumes,    and   so    it   seems    as    if  a 
wire  should  produce  more   current.      But   the  problem   is   complicated   by   the 
electric  field  distribution  which  determines    the    longitudinal    dimension   and 
a     the  multiplication   coefficient.      As  will    be  seen    total    current    is   not 
a    function  of     a     and     6      alone,    but   a   function  of   the  product      (ao)      or 
(a*6*)    in   the  flow  situation. 

2.      Total    Current   Considerations 

As   pointed  out    in   Part    II,    Section   3,    since   the   total    current    at 
the  anode   surface    is   due  to  electrons  entering   the  anode   from   the 
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multiplication  region  (the  positive  ions  being  repelled),  the  total  discharge 
current  can  be  expressed  as  a  function  of  the  number  of  electrons  entering 
and  their  drift  velocity  (the  assumption  is  made  that  if  negative  ions  are 
present,  they  give  up  their  electron,  soon  after  drifting  into  the  high 
field  region;  this  assumption  is  supported  by  Reference  15  and  some 
sample  calculations  presented  in  Reference  8). 

A  single  electron  entering  the  high  field  region  will  exponentially 
multiply  to  produce 

J      adx 
x 
n  =  n  e 
o 

electrons  at  the  anode  surface  with  n  =  1  .   Since  the  current  density  is 

o  T 

the  current  per  unit  area, 

j  =  en  v  =  i/A  (50) 

Je     e  e 

the  total  current  (which  must  be  constant  across  the  discharge)  can  be 

expressed  as 

f     adx 

I  -  '>       =  jeA  (51) 

This  allows   the   total    current   to  be  discussed    in    terms   of  electrons 

arriving   at    the  anode.      That  number   is    strongly   dependent    upon   how  many 

electrons  enter   the  multiplication    region,  whether  as    free  electrons,   or 

attached  to   atoms  or  molecules   as   negative    ions.      The  electron    input    is 

described  as      i         in    the   above  equation   and   is 
o  ^ 

i      =(en      v     +env)AH  v-e  (52) 

o  e     e  -   - 

o 


kO 


where 


n     =  free  electron  density  from  drift  region 
o 

n_   =  negative  ion  density  from  drift  region 
v  ,v.  =  particle  drift  velocities 

A    E  area  through  which  particles  pass  into  multiplication  region 
Y    =  secondary  electron  supply  rate. 

Without  covering  many  details  or  methods  already  discussed  in 
Section  A,  the  aoove  considerations  lead  to  the  normalized  current  equation 
below 

i       2Ca  6  6"n(a*<5*-l) 
i*  -  j--  y,ce  (5  3) 

o 


where 


v*  z    i   /i    -   normalized  secondary  electron  supply  rate 
'     ou  oo  7  KK  ' 


From  equation  52  it  can  be  seen  that  the  units  of  y  are   electrons/second, 
and  that  it  is  an  electron  flow  rate  into  the  multiplication  region.   It  is 
the  counterpart  of  the  Townsend  secondary  ionization  coefficient  in  a  non- 
flowing  discharge. 

Understanding  the  nature  of  y  ,   we  are   nrow  able  to  analyze 
possible  flow  effects.   Electrons  coming  into  the  ionization  region  can 
arrive  in  two  forms,  either  as  free  electrons  at  the  characteristic 
electron  drift  velocity  which  is  generally  much  higher  than  typical 
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flow  velocities,  or  attached  to  a  heavier  particle  in  the  form  of  negative 
ions,  etc.,  whose  drift  velocities  and  therefore  delivery  rates  can  be 
greatly  altered  by  the  flow. 

Now  we  discuss  the  third  component  of  current,  the  positive  ion 
contribution.   In  the  determination  of  total  current  by  counting  electrons 
flying  past  the  anode  surface,  the  positive  ions  play  an  indirect  role  by 
supplying  secondary  electrons  primarily  through  impact  with  the  cathode  where 
energetic  col  lions  can  result  not  only  in  neutralization  of  the  ions,  but 
also  release  of  additional  electrons  which  then  begin  to  move  toward  the 
anode.  These  secondary  electrons  may  proceed  as  free  electrons  or  attach 
to  particles  and  proceed  as  negative  ions,  thus  the  positive  ion  contri- 
bution shows  up  as  an  enhancement  of  both  y   and  y   •   Since  both 
r  'e       'n- 

positive  and  negative  ion  drift  velocities  are  of  the  same  order  of  magnitude 
as  the  flows  being  discussed,  it  can  be  anticipated  that  the  total  current 
would  be  significantly  affected  by  flow  through  ion  slip  -  and  electron 
trapping  discussed  in  Part  I IA. 

Considering  the  free  electron  component  alone,  the  total  input  of 
electrons  to  the  multiplication  zone  is,  from  equation  52 


i"  =  (e  n  v  )A  (54) 

o      e  e 


At  first  it  would  seem  that  electron  velocity  in  the  drift  region  should 
not  be  significantly  affected  by  the  flow  due  to  low  electron  slip  and 
negligible  trapping  due  to  space  charges.   However 


ve  s  keE(x)  (55) 


kl 


and 


e  k     E 

u  e      u 


V-T^-irr1-6*  (56) 

e  e     o 

o 


In   nearly  all    non-uniform  geometries,    the  electric   field    is   given   by  an 
expression  of    the   following    form  - 


E(*>    -r^irL^r\  (57) 


x    log(CONST) 


with    the    result   that 


E*(x)    =  V' 


and  at  breakdown, 


Eb*  =  Vb*  (58) 

Therefore   from  equations   5^,   55,    and  58  - 


,*.^ou     =  Veu  =  Vb-,n    *  (59) 

o        i-  n      v  e 

oo  eo  eo 


This  equation  shows  that  even  free  electron  velocities  can  be  altered 
indirectly  by  flow  modifications  which  affect  the  breakdown  voltage 
and  hence  the  electric  field  even  in  regions  where  electron  trapping  is 
not  a  significant  factor.    Vb"  and  therefore  electron  speed  can  be  altered 
by  30%  or  more  at  moderate  flow  velocities. 

Now  the  total  current  can  be  written  as  follows  - 


^3 


Ca  <5  5*n(a*<5*-l) 


i*  =  n  *Vb*e     °  °  (60) 


The     n   *     term   is  a   function  of   the   total    current,   and  all    associated 
e 

secondary  electron   supply  processes   such  as   cathode  emission,    photoioniza- 

tion,  etc.,  whose  variations  with  flow  are  as  yet   undetermined. 

In   a  discharge  where  many  of  the  secondary  electrons   are   liberated 

by  positive    ion    impact  with    the  cathode,   however,     n    *     will    be  a   function 

e 

of  the  ion  slip  and  the  current  equation  may  assume  the  following  form  - 


i*,Vb*  0+S)^Y*e  °  °  (61) 


Neglecting     y*     temporarily,   Equation  61   will    be  used   to  qualita- 
tively examine   the  behavior  of   the  breakdown   current  with   flow  and   turbu- 
lence and   to  attempt    to  explain   some  of   the  observed  differences    in    the 
point-to-plane  and  wire-to-plane     Vb*     and      lb     vs.      velocity   data  collec- 
ted during   the   current    research  effort. 

Figures  16,  17  and  18  are  typical  plots  of  breakdown  voltage  and 
current  vs.  velocity  for  the  point  and  wire  configurations  and  will  be 
referenced   frequently    in   the   following   discussion. 

a.      Point-Plane   Current    Interpretations 

Since  breakdown   voltage  variations  with   flow  were  already 
discussed   in   detail    in   Part    I  I. A. 5,  we  will    concentrate  on   some   typical 
point-plane   current  data   first.      Beginning  at    the  origin      (U=0)    in   Fig. 
17a,    it   can   be   seen   that,   as   velocity    increases    in   the  normal    flow 
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(positive)    direction,    the   current  begins    to   rise   (region    1),    reaches   a    low 
maximum   (region   2)    and   then   declines   slightly    (region    3).      A   typical 
Vb"      curve    (Fig.    8)    shows  a   depressed  breakdown  voltage    in    this    region, 
which   from  Equation  61    tends    to  decrease      lb      but    the   increasing   velocity 
term  helps    to   compensate  for   the   reduction   as   does   the  product      (a*5*)      in 
the  exponential,   and   the   current    in    regions    1    thru   3    is   dominated  by   the 
behavior  of     5*      (in    this   velocity    region    turbulence    intensity    is   not   yet 
sufficient    to  cause  significant  enhancement  of     a*  ).      As   velocity    increases, 
if  an  effective   turbulence  generator   is    in   place,      a*     becomes   dominant    due 
to   its    rapid  exponential   growth  with    increasing   velocity    (Equation   kQ) ,    and 
in   accordance  with  Equation  61    an  overexponent i al    rise      in      lb      is   noted 
(region   k  of   Fig.    17a).      The   figure   also   shows     Vb*      leveling   at    the  normal 
flow  asymptote  which  would   also  cause  a    leveling  of   the      lb      characteristic 
or  at   least   a    leaning    toward   a      (1+S)      linear   rise    in   accordance  with   Equa- 
tion  61. 

b.      Reverse   Flow  Current 

Interpretation  of    the   reverse  flow   (negative   velocity)    current 
and  voltage   characteristics    in    the  point-plane   and  wire-plane    is   essentially 
the   same,    so   the  following  discussion   applies   to   both   configurations.      Since 
no   turbulence  was    introduced    in    the   reverse  flow  mode,   other   processes    tend 
to   dominate  breakdown   characteristics    in    this    regime.      Equation   kS   shows   a 
rapid  decrease    in     5*     with    reverse  flow  which   causes   an   equally    rapid    rise 
in      Vb*      ;    the  Vb*      increase  and      (1+S)      velocity   term  decrease  might   be 
expected  to  neutralize  each  other  as    in   the  normal    flow  mode,    leaving    the 
major    influence  on      lb      the     5*      term  again.      Since     5*      is    rapidly   decreas- 
ing,   this    line  of    reasoning  would  have      lb      decreasing   at   a   fairly    rapid 


rate.      However,    the  experimental    data  consistently   show  a   steep   rise 
in      lb      at   relatively    low  negative  velocities    (region   5  of   Fig.    17a),    and 
a    leveling   or  decrease  at   higher  velocities    (regions   6   and   7).      The  high 
velocity   reduction    in      lb      is  easily   understood  via   Equation   61,    since    the 
velocity   term      (1+S)      will    eventually  approach   zero  at   high  enough  nega- 
tive velocities   causing      lb      to  approach   zero.      Also  since     6*     is    con- 
stantly decreasing  with    increasing    reverse  flow,    this    tends    to   reduce      lb 
more  at  high  velocity.      (In  one  experimental    case      lb     was    reduced    to 
approximately      1/10     of   its  no   flow  value).      It   seems   that   the  only   factor 
going    in   favor  of   increasing   reverse  flow  current    in    the   low  velocity 
region    is   the   increase   in      Vb*   ;      but  numerical    calculations    show   that    the 
velocity    term  more   than  offsets   this    factor.      Careful    review  of   the  equa- 
tions  brings   aid   from  an    unexpected  source   to  explain    this   strong    low 
negative  velocity   current  enhancement   through   the     a*     term. 

Recalling    the  form  of     a     from  Equation   6,    it   can   be  seen   that 
a      is   a   function  of     E      as  well    as      P    ,    and  as      Vb*      increases    in    the 
reverse   flow  mode,      E      also   increases,    since   from  Equation    1      E  -   V   . 
Mathmati  cal  1y   this    leads   to    the   following  enhancement  of     a     due   to   the 
rise  of     E      for   the  point-plane  - 

ct*(AE)    z  exp{(-BP/Vbojf(R   log(  WR)](1/Vb*-1)  }  (62) 

Another   subtle  effect  may  also  come   into  play  here,   since,   as   flow   is 
directed   into  the  quasi-steady  no-flow  anode  extension    region,    the   initial 
hydrodynami c  effect  will    be   to  promote  a   flattening   and   radial    spreading 
of  the  dense  space   charge  boundary.      In    this   case   the  point    radius      R 
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should  be   replaced  by      R__       which  will    be  a   complicated   function   of   the 
flow  and  electric  fields.      Qualitatively,   however,    it    is  expected   that   after 
initial    low  velocity  spreading  which  enhances     a*      (via   Equation  62)    and 
therefore      i*   ,      i*     will   eventually  decline   due    to  the  continued  decreases 
in     5*     and    the      (1+S)      velocity   term   in   accordance  with   Equation  61    and 
R  _P     will    also  decline.      The  end   result  of   these  processes    is   a    rapid 
reverse   flow  current   rise  and   fall    as    illustrated   by   Fig.    Ja.        These   two 
reverse  flow  effects        a*(AE)      and     R  •*-   R___        are  not   peculiar   to    reverse 
flow,   but    in    the  normal   mode   they  must    compete  with  other   strong    influences 
which  may   tend    to  mask   their  effects  except    under  special    circumstances. 
For  example,    in    the  normal   mode,   gas   flow  and  positive    ion    repulsion   from 
the  anode   act    in    the   same   direction  on    the  high   density   positive   space 
charge,   so  gas   convection   by    itself  does  not   promote    radial    spreading   as 
it   does    in    the   reverse   flow  mode  where  electrostatic   repulsion   and   convec- 
tion oppose  each  other.      If,    however,    the  normal    flow  becomes   highly   tur- 
bulent,   the    random  gas   velocity   and   density   fluctuations    in    the    ionization 
region    tend   to  enhance    ionization   and   the  transport  properties   of    the  gas 
resulting    in   faster   radial    diffusion    and  dynamic  spreading   of   the  positive 
charges,    and   as    in    the    reverse   flow  case,    the  effective    radius   of   the 
virtual    anode    is    increased.      Reference   3  presented   photographs   of  various 
discharge   conditions  which  visually    illustrates   this    turbulent   enhancement, 
the  glow   region    in    the   turbulent    case  being   3*5   times   the   radius  of   the 
no- flow  or    laminar   flow  cases. 

c.      Wire-Plane  Normal    Flow  Breakdown   Current   and  Voltage 

Breakdown  data  for   this   configuration  was    limited    to  a   short 
gap    length    range  above   the  existence  point   for  the   corona   discharge      (~2   cm) 
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due   to  power   supply   problems  which    limited  available   voltage   to   approxi- 
mately     33KV.      This    limitation   effectively    restricted  wire-plane   gap    length 
to     2.8  cm     or    less  with    flow.      Since    this    is    less    than      1    cm      into   the 
corona  existence    region,    it    is   not  possible   to  give   a   full    accounting  of 
flow  effects.      In   general,    however,   non-turbulent   curves  were   usually 
flat  with    little  enhancement  of   current,   while  most    curves  with    turbulence 
screen   #9    in   place   showed  slow  but   steady    improvement  with   velocity.      One 
interesting    test  was   performed  with   a  gap    length   below   the   corona  existence 
point,    and   the    introduction   of    turbulent   flow   resulted    in    the   flow  of   a 
small    current  before   breakdown.      This    is  most    likely   due   to  the   pressure/ 
density   fluctuations    in    the   vicinity  of    the  wires  which   cause    corresponding 
fluctuations    in      E/P    .      If  the    lowest   pressure   variations    result    in    a  high 
enough  value  of    the    reduced  field      (E/P),    some    ionization  will    begin    to 
occur  before   breakdown. 

The  wire-plane   breakdown   voltage   curves  were  mentioned    in   Section    I 
A. 6   and  the   general    trends   shown    in    Figs.    8  and    18.      One   additional    obser- 
vation will    be  made  here.       It  has   been   noted    in    the   point-plane   configura- 
tion   in    the  normal    flow  mode,    that  as   gap    length    is   shortened   toward    the 
corona  existence   point,    breakdown   voltage  will    rise   at    lower   flow  veloci- 
ties.     In    the  most  extreme  point-plane   case,      Vb*     had   reached   the  normal 
flow  asymptote   at     kQ  m/s     with  no_  turbulence.      In    fact    the   curve    looked 
almost    like   a  wire-plane   curve.      This   seems    to    infer   that   if  gap    length    is 
in    the   corona  existence    region,    but   close   to   the   existence   point,    the 
normal    flow     Vb*     characteristic  may   tend   toward   the  wire-plane   shape  of 
Fig.    18.      That    is,    if   the  gap    is   short    relative   to   the   corona  existence 
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point  a  low  velocity  rise  in  Vb  will  occur.   Additional  data  at  various 
gap  lengths  are  necessary  to  confirm  this  preliminary  observation.   If, 
in  fact,  this  is  a  true  characteristic,  it  is  most  likely  related  to  the 
variation  of  the  axial  breakdown  electric  field  distribution  as  the  gap 
length  is  reduced. 

The  remaining  sections  review  requirements  for  the  optimum 
pumping  of  a  laser  gas  mixture  vs.  electrode  geometry,  a  final  over- 
view of  data  and  theory  developed  earlier,  and  some  conclusions  and 
recommendations  for  further  research. 
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C.      Electrode   Design/Optimization 

In  electric  discharges   designed   for    laser  pumping,    two   basic   require- 
ments must   be  satisfiedo      First,    the  electrodes  must  efficiently   produce 
a   sufficient  electron   population    to  permit  effective  energy   transfer   to 
the    laser  gas.      Secondly,    the  electrode   design  must   be   such    that    the 
reduced  electric  field    (E/N  or  E/P)    in    the  pumping    region   optimizes  elec- 
tron energy   for  maximum  energy    input    to   the   desired  vibrational    or 
electronic  energy    levels*      Degradation    in   either  of    these  critical    areas 
will    result    in    reduced    laser   power  output. 

Discharge   cavity  size    limitations   or  specifications  will    most    likely 
dictate   the  maximum  gap   spacing  between   anode  and   cathode,    so  within    that 
limitation    the   combination  of    the   stressed  electrode    radius   and  gap    length 
must   be   determined  which   best   satisfies   the   requirements   above.      Since    the 
E-field   distribution    is   a   function  of  both   parameters,    they   cannot    be 
considered   separately.      In   addition,    it   has   been   shown   that    space   charge 
effects   via   the  "anode  extension"   and   radial    spreading   cause   significant 
changes    in    the  "apparent"  electrode  geometry  which   should   be   taken    into 
account    in    the   final    design.      The  gasdynami c   interactions  with    these   space 
charge    induced  dimensional    changes,    in   effect,    provides   an   extra   degree   of 
freedom   in   a   fixed  geometry   cavity   through   some   control   of   the    radial 
spreading  and  anode  extension.      This    in    turn   permits,    to   some  extent,    ad- 
justment of     E/N      and  electron   production   via   flow  and   turbulence. 
Additional    benefits   arise   through   the  enhancement  of   transport   properties, 
improving   stability  while  concurrently    increasing  electron   production. 

Calculation   of   the    reduced  fields   via   the  adaptation   of   the  hyper- 
boloidal    point-plane,   and  an   assumed  wire-plane   distribution    indicates 


E/N   values  in  the  drift  regions  of  the  discharges  tested  slightly  higher 
than  the  optimum  range  specified  in  Reference  1,  a  condition  which  can  be 
remedied  by  proper  adjustment  of  the  geometric  parameters,  provided  that 
an  adequate  power  supply  is  available. 
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III.      SUMMARY 

A.        Data 

The     Vb-     data   for  both   configurations   are  quite   similar   in    the 
reverse   flow  modes,    showing   a   fast   rise   to  an   apparent   asymptotic    limit 
(Figs.    16   a,b).      This    is   presumed  due   to  collapse  of   the  anode  extension 
by  the  gas   flow. 

The  normal    flow  mode   shows     Vb*      initially   decreasing    for  all    gaps 
and  geometries  with   an  eventual    upswing  which   depends   strongly  on    turbu- 
lence   intensity,   and   appears    to   also  depend  on   gap    length    relative   to   the 
corona  existence   point,   with  early    (low  velocity)    rise   strongly    favored    in 
the  high   turbulence   and/or   short   gap   cases    (Figs.    8,    16   a,    b,    18). 

The  breakdown   current  data    in    the   reverse   flow   illustrated   similar 
behavior    in    the  point  and  wire-plane   configurations,    starting  with   a   sub- 
stantia]   low  velocity   rise   followed  by  a    decline   at  higher   velocities. 

The  point-plane  normal    flow  current   behavior  exhibits    several    features 
under  various    conditions.      The  most   promient  of   these  are  outlined    in 
Fig.    17a.      Many   smaller  scale  perturbations   are  evident    in    some    individual 
sets  of  data,    but    the  high    turbulence,    higher   current    curves   seem   to 
follow   the   same  basic  pattern    in    the  point-plane   configuration.      The  wire- 
plane   current  curves  exhibit  much  more   scatter  and   trends   are   not   as  well 
defined  as  with   the  point  geometry.      This    is   probably  attributable   to    the 
fact   that   data  were  not   accessible    in   a  wide    range  of  gap    lengths,    but 
were   restricted   to  gaps  close   to   corona    initiation  where      Vb     and     Vc     are 
not   far  apart,   and   therefore   current   levels   are    relatively    low. 
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B.        Theory 

The    theory   developed   in   Part    II    accounts   for  the  major   features   of 
the     Vb     and      lb      curves   described  above  as   can   be  seen    in    the  various 
figures,   but  many  of   the  analytical    terms   can   be   improved  and    it    is 
believed   that   further  study  will    lead   to  analytic  expressions    for   the 
various   parameters  of   the  theory  which   have   been  mentioned.      The  clear 
physical    basis  of   the   theory    lends   to  easy    interpretation  of  various    flow 
related  phenomena   and   should  prove   useful    in   visualization   and   design   of 
flow-type  systems   and  possibly    in   trouble-shooting  of  existing   flow   system 
p  rob  1  ems . 
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IV,   CONCLUSIONS  AND  RECOMMENDATIONS 

A.  Conclusions 

1.  The  basic  anode  extension  model    developed  herein   appears    to 
offer   a   simple  approach   to   the    interpretation   and   calculation   of  basic 
gasdynamic  effects  on   high  pressure   corona  discharges. 

2.  The   flow  field   interactions   discussed  are   basic,    so   the  approach 
used    in  modelling    the   positive  point-plane  electrode   system  should  be 
easily  applied   to  other  gap   configurations,    polarities,    gas   mixtures,    and 
pressure   regimes. 

3.  A  more    rigorous  evaluation   of  various    terms    in    the  basic  model 
should  provide  better   insight    into   the   dominant   forces    in   various    flow 
regimes,   and  most    likely  will    result    in   definition   of   some  of   the  small 
scale   features  not   yet  well   explained. 

B.  Recommendations 

1.  Continue  data   collection   on    the  wire-plane  electrode   configura- 
tion when   a    larger  power   supply    is   available    (50KV,   50ma   should   be 
sufficient  with   present   hardware). 

2.  Develop  a   good  analytical   expression    for   the  effective   electrode 
radius   for    insertion    into   the  present   point-plane   computational   model. 

3.  Include     a*(AE)    and    radial    spreading  effects    in   normal    flow 
computational    scheme. 

k.      Test   point-plane    in    reverse   flow  mode  with   turbulence   at    low 
velocities    to  verify  Equation  45. 

5.      Develop  comprehensive  analysis  of   flow  on     y*     for  various    types 
of  gas  mixtures,    i.e.,   attaching,   non-attaching,   etc. 
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GRID 


(a  )      Normal  Flow 


£> 


CATHODE 


ANODE 


(b)      Reverse  Flow 


GRID 

CATHODE 

$>\ 

•                                                  ANODE 

(c)     Cross -Flow 

Figure  2.    VARIOUS   ORIENTATIONS  OF  ELECTRODES 
AND   TURBULENCE  GRIDS  TESTED. 
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opp 
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Figure  4.  ELECTRIC  FIELD  AND  VOLTAGE  DISTRI 
BUTION  BETWEEN  PARALLEL  PLATES. 
Subscripts1    ( I  )  No  Space  Charge 

(2)  Moderate  Space  Charge 

(3)  High  Space  Charge 
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ANODE  PINS 


ROW    1 
(a) 


Equipotentials 


ROW   2 
(b) 


-XrJ\  - 


CATHODE 


ROW    3 
(c) 


1  .    Applied  Field 

2.  Space  Charge  Field 

3.  Total  Field 

Figure   5.    (a)  Schematic  of   Field  Lines  and  Equipotentials. 

(b)  Apparent  Dimensional  Changes  Due  to  Space 
Charge. 

(c)  Applied,  Space  Charge  and  Total  Field 
Distributions. 

X     and  x'    are  Critical  Distances  From 
c  c 

Anode   Where  Electron  Avalanching  Begins 

(E  =  E.  )  With  and  Without   Space  Charge 
i 

Effects. 
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SECTION 


CATHODE 


SECTION  A-A 


(b) 

Figure  10.      (a)    Activity  in  a  Parallel  Plate  Gap. 

(b)    Modified  E- Field  Due  to  Space  Charge. 
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